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Abstract 

Boltzmann equations are often used to study the thermal evolution of particle reaction networks. 
Prominent examples are the computation of the baryon asymmetry of the universe and the evolution 
of the quark-gluon plasma after relativistic heavy ion collisions. However, Boltzmann equations 
are only a classical approximation of the quantum thermalization process which is described by 
the so-called Kadanoff-Baym equations. This raises the question how reliable Boltzmann equations 
are as approximations to the full Kadanoff-Baym equations. Therefore, we present in this paper a 
detailed comparison between the Kadanoff-Baym and Boltzmann equations in the framework of a 
scalar d> 4 quantum field theory in 3+1 space-time dimensions. The obtained numerical solutions 
reveal significant discrepancies in the results predicted by both types of equations. Apart from 
quantitative discrepancies, on a qualitative level the universality respected by the Kadanoff-Baym 
equations is severely restricted in the case of Boltzmann equations. Furthermore, the Kadanoff- 
Baym equations strongly separate the time scales between kinetic and chemical equilibration. This 
separation of time scales is absent for the Boltzmann equation. 
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I. INTRODUCTION 


One of the most attractive frameworks to explain the matter-antimatter asymmetry of the 
universe is the so-called leptogenesis mechanism [1-3]. Here, lepton number violating inter¬ 
actions in the early universe produce a lepton asymmetry which is subsequently converted 
to the observed baryon asymmetry by so-called sphalcron processes. For the dynamical 
generation of the lepton asymmetry it is necessary, that the universe was in a state out 
of thermal equilibrium [4], The standard means to deal with this nonequilibrium situation 
are Boltzmann equations. However, it is well known that (classical) Boltzmann equations 
suffer from several shortcomings as compared to their quantum mechanical generalizations, 
the so-called Kadanoff-Baym equations. This motivates a comparison of Boltzmann and 
Kadanoff-Baym equations in order to assess the reliability of quantitative predictions of 
leptogenesis scenarios. 

In addition to leptogenesis, there are various other systems which warrant a comparison 
between Boltzmann and Kadanoff-Baym equations: In particular, a strong motivation is 
furnished by relativistic heavy ion collision experiments which aim at testing the quark- 
gluon plasma. In these experiments the quark-gluon plasma is produced in a state far from 
equilibrium. Recently, however, experiments claimed that the approach to thermal equi¬ 
librium should happen very fast, and that the evolution of the quark-gluon plasma could 
even be described by hydrodynamic equations [5-8], which arise as approximations to Boltz¬ 
mann equations. In this context it is important to note that different quantities effectively 
thermalize on different time scales [9]. Thus, one might face the situation that, although 
the full approach to thermal equilibrium takes a very long time, certain quantities, which 
are sufficient to describe the quark-gluon plasma with hydrodynamic equations, approach 
equilibrium values on much shorter time scales. 

In order to derive Boltzmann equations from Kadanoff-Baym equations 1 , one has to em¬ 
ploy several approximations, among them a first-order gradient expansion, a Wigner trans¬ 
formation and a quasi-particle (or on-shell) approximation [12-16]. However, it is known, 
that the gradient expansion cannot be justified for early times. Consequently, one might 
expect that Boltzmann equations fail to describe the early-time evolution and that errors 
accumulated for early times cannot be remedied at late times. Of course, a Wigner trans¬ 
formation itself is not at all an approximation, but in order to make it available, one has to 
send the initial time to the remote past. Whereas Boltzmann equations imply the assump¬ 
tion of molecular chaos, meaning that two particles were uncorrelated before their collision, 
Kadanoff-Baym equations take these memory effects into account. Numerical solutions of 
Kadanoff-Baym equations revealed that this memory is lost gradually. Consequently, for 
late times it is indeed justifiable to send the initial time to the remote past. However, for 
early times this is certainly not the case. Furthermore, as a consequence of the quasi-particle 
approximation, the conservation of momentum and energy prevents Boltzmann equations 
from describing thermalization in 1 + 1 space-time dimensions. In contrast to this, it has 
been shown in the framework of a scalar <f> 4 quantum held theory that this is feasible with 
Kadanoff-Baym equations [17]. The reason for this qualitative discrepancy is that Kadanoff- 
Baym equations take off-shell effects into account [18], which are neglected in Boltzmann 
equations. Of course, in 3+1 dimensions both types of equations are capable of describing 
thermalization. In the case of leptogenesis, however, the on-shell character of the Boltzmann 


The connection between Boltzmann equations and classical field theory has been treated in Refs. [10, 11]. 
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equation leads to a further inconsistency: All leptogenesis scenarios share the fact that some 
heavy particles decay out of thermal equilibrium into the particles which we observe in the 
universe today. The spectral function of a particle that can decay into other particles is 
given by a Breit-Wigner curve with a non-vanishing width. By employing the quasi-particle 
approximation we reduce this decay width of the particles to zero, i.e. a Boltzmann equation 
can only describe systems consisting of stable, or at least very long-lived, particles. After all, 
how does the on-shell character of the Boltzmann equation affect the description of quantum 
fields out of thermal equilibrium in 3 + 1 dimensions? 

When applying Boltzmann equations to the description of leptogenesis, the standard 
technique to construct the collision integrals — before employing further approximations 
- is to take the usual bosonic and fermionic statistical gain and loss terms multiplied 
with the S-matrix element for the respective reaction [19, 20]. These S-matrix elements are 
computed in vacuum, and one may wonder of which significance they are for a statistical 
quantum mechanical system. 

All these shortcomings of Boltzmann equations lead to the conclusion that one should 
perform a detailed comparison between Boltzmann and Kadanoff-Baym equations [13, 21- 
23], such that one can explicitely see how large the quantum mechanical corrections are. Due 
to the complexity of the problem, we restrict ourselves for the moment to a $ 4 quantum field 
theory in 3 +1 space-time dimensions. Of course, in this framework one can neither describe 
the phenomenon of leptogenesis nor thermalization after a heavy ion collision. Nevertheless, 
it may well serve as starting point for further research, and certainly permits to present a 
detailed comparison between Boltzmann and Kadanoff-Baym equations, which may point 
to interesting phenomena to be investigated in more realistic theories. 

In general, when studying systems out of thermal equilibrium by means of Kadanoff-Baym 
equations, it is crucial to start from a <h-derivable approximation, since these approximations 
ensure the conservation of energy and global charges [24-26]. The 2PI effective action 
furnishes such a ^-derivable approximation [27-29] and has proven to be an efficient and 
reliable tool for the description of out-of-equilibrium quantum fields in numerous previous 
treatments [17, 30-34]. In this work, we start from the 2PI effective action truncated at 
three-loop order. The Kadanoff-Baym equations can be obtained by requiring that the 
2PI effective action be stationary with respect to variations of the full connected two-point 
function [17, 18]. In order to derive the corresponding Boltzmann equation, subsequently one 
has to employ a gradient expansion, a Wigner transformation, the Kadanoff-Baym ansatz 
and the quasi-particle approximation [12-16]. While the Boltzmann equation describes the 
time evolution of the particle number distribution, the Kadanoff-Baym equations describe 
the evolution of the full quantum mechanical two-point function of the system. However, one 
can define an effective particle number distribution which is given by the full propagator and 
its time derivatives evaluated at equal times [17, 31]. Finally, we solve the Boltzmann and 
the Kadanoff-Baym equations numerically for spatially homogeneous and isotropic systems 
in 3+1 dimensions and compare their predictions on the evolution of these systems for 
various initial conditions. 
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FIG. 1: Closed real-time path C. This time path was invented by Schwinger [35] (see also [36, 37]) 
and applied to nonequilibrium problems by Keldysh [38]. In order to avoid the doubling of the 
degrees of freedom, we use the form presented in Ref. [13]. 


II. 2PI EFFECTIVE ACTION 


In this work we consider a real scalar quantum field, whose dynamics is determined by 
the Lagrangian density 

se = -l(a„4>) (a"t) - - It 4 . 

The minus sign of the kinetic term indicates that we use the metric where the time-time 
component is negative. As we will compute the evolution of the two-point Green’s function 
for a nonequilibrium many body system, already the classical action has to be defined on the 
closed Schwinger-Keldysh real-time contour, shown in Fig. 1. The free inverse propagator 
can then be read off the free part of the classical action 


Io = ~ 


1 

2 


J d 4 x j d A y [<f> (x) G 0 1 (x, y ) <h (y)] 
c c 


> 


where 

Gq 1 ( x , y) = (d xl *d yti + m 2 B ) 5 C (x - y) . (1) 

We consider a system without symmetry breaking, i.e. ($ (x)) = 0. In this case the full 
connected Schwinger-Keldysh propagator is given by 


G (x, y) — (Tc (x) $ (y)}) . 

Accordingly, for Gaussian initial conditions the 2PI effective action can be parameterized in 
the form [17, 27-29] 

r [G] = ^tr c log c [dA 1 ] - itr c [Gq 1 G] + T 2 [G] + const . 

iY 2 [G] is the sum of all two-particle irreducible vacuum diagrams, where internal lines 
represent the full connected propagator G(x,y). Of course, for an interacting theory we 
cannot compute T 2 [G] completely, and we have to rely on approximations. In this work we 
apply the loop expansion of the 2PI effective action up to three-loop order. The diagrams 
contributing to T 2 [G] in this approximation are shown in Fig. 2. We find [18]: 

r 2 [G] = —— J d 4 x[G (x,x) G (x,x)] 
c 

+ ^ j d A x j d 4 y [G (x, y) G (x, y) G (y, x) G (y, x)\ . 
c c 
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FIG. 2: Two- and three-loop contribution to T 2 [G]. The lines represent the full connected 
Schwinger-Keldysh propagator. 

III. KADANOFF-BAYM EQUATIONS 

The equation of motion for the full propagator reads [27, 28] 

jr [g] 

SG(y,x) 

It is equivalent to the Schwinger-Dyson equation 

G- 1 (x, y ) = iGq 1 ( x , y)-U (x, y) , (2) 

where the proper self energy is given by 

n (x, y) = 2i ^ 2 ^ = -^Sc (x — y)G (x, x) - ~-G (x, y) G (x, y) G (y, x) . (3) 

dG[y,x) 2 o 

After we have inserted the free inverse propagator given in Eq. (1), we convolute the 
Schwinger-Dyson equation (2) with G from the right: 

i ( d x pd x ^ + m%) G (x, y) = 5 C (x - y) + J d A z [II (x, z ) G (z, y)] (4) 

c 


Next, we define the spectral function 2 

G g (x,y) =i([$ (V) ,^(y)]_> 

and the statistical propagator 3 

Gf (x,y) = ^ ([$ (x) ,$(y)] + ) 
such that we can write the full propagator as 

G (x, y) = G f (x, y) - ^sign c (x° - y°) G e (x, y) . (5) 


2 From the definition of the spectral function we see that it is antisymmetric in the sense that G e (x , y) = 

—G e (y,x). Furthermore, the canonical equal-time commutation relations give {G e (x,y)) x0 0 = 0 and 

(d y oG e (x, y)) x o =y o = -$ 3 O - y)- 

3 In contrast to the spectral function, the statistical propagator is symmetric in the sense that Gf (x, y) = 
Gf (y,x). 


5 




FIG. 3: One- and two-loop contribution to the proper self-energy II. Again, internal lines represent 
the full connected Schwinger-Keldysh propagator. The tadpole represents the local part which 
causes a mass shift only. The setting-sun diagram represents the nonlocal part and leads to 
thermalization. 


Note that for real scalar quantum fields both the statistical propagator and the spectral 
function are real-valued functions [17]. The spectral function describes the particle spectrum 
of our theory. From its Wigner transform we can obtain the thermal mass and the decay 
width of the particles in our system. On the other hand we will define an effective particle 
number density given by the statistical propagator and its time derivatives evaluated at 
equal times. From Eq. (3) (and Fig. 3) we see that the self energy contains a local and a 
nonlocal part: 

n (x, y ) = -iS c (x - y ) Ii (local) ( x) + Yi^iocai) (x, y) . 

The local part of the self energy only causes a mass shift, which can be included in an 

effective mass: 

M 2 (x) = m B + U {local) ( x) = m 2 B + ^G F (x, x) . (6) 

After inserting Eq. (5) into Eq. (3), we can decompose the nonlocal part of the self energy 
in exactly the same way as we did for the propagator: 

n ( nonlocal ) ^ = Up _ l s i g n c (a; 0 - y°) U e (. X, y) . 

We find 

A 2 / 3 \ 

n F (x, y) = - — \G f (x, y) G F (x, y) G F (x, y) - -G e (x, y) G e (x, y) G F (x, y) J 

and 

n e (x, y) = —— (3GV (x, y) G F (x, y) G e (x, y) - -G e (x, y) G e {x, y) G Q {x, y) J . 

When we insert all these definitions into Eq. (4), we observe that it splits into two com¬ 

plementary real-valued evolution equations for the statistical propagator and the spectral 
function, respectively [17]. These are the so-called Kadanoff-Baym equations: 

y° x° 

{-d x nd x ^ + M 2 (a;)) G F (x, y) = j d A z If F (x, z) G e {z, y) - J d A z U e (x, z) G F (z, y) (7) 

o o 

and 

x° 

{-d x »d Xii + M 2 (x)) G 8 (x, y) = - J d A z U g (x, z) G Q (z, y) . (8) 
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For a spatially homogeneous system, one can Fourier transform these equations with respect 
to the spatial relative coordinate. Furthermore, in an isotropic system the propagator will 
depend only on the modulus of the momentum. As explained in more detail in Refs. [17, 31], 
one can define effective kinetic energy and particle number densities uj ( t , p) and n ( t , p) which 
are given by 


uj 2 ( t,p) 


( d x od y oG F (x°,y°,p) \ 
v G F (x 0 ,y 0 lP ) ) x o =y0=t 


(9) 


and 


n(t,p) = u (t,p) G f ( t,t,p) 


1 

2 ' 


( 10 ) 


However, we stress that the Kadanoff-Baym equations are self-consistent evolution equations 
for the full propagator of our system, and that one has to follow the evolution of the two- 
point function throughout the whole x°-y°-plane (of course, constrained to the part with 
a: 0 > 0 and y° > 0). One can then follow the evolution of the effective particle number 
density along the bisecting line of this plane. 

We would like to emphasize that the only approximation involved in the numerical solu¬ 
tion of the Kadanoff-Baym equations is the loop expansion of the 2PI effective action. In the 
next section we will describe the approximations which are necessary to derive a Boltzmann 
equation from the Kadanoff-Baym equation for the statistical propagator (7). 


IV. BOLTZMANN EQUATIONS 

It is well known how Boltzmann equations can be obtained as an approximation of the 
Kadanoff-Baym equations [12, 13, 16]. In this section we briefly review the standard deriva¬ 
tion: One has to employ a Wigner transformation, a gradient expansion, the Kadanoff-Baym 
ansatz and the quasi-particle approximation. 

First, we subtract the Hermitian adjoint of Eq. (7) from Eq. (7) and re-parameterize the 
propagator and the self energy by center and relative coordinates 

G (u, v) = G y—, u-vj . 

Next, we define X = and s = x — y, and observe on the left hand side of the difference 
equation that 

d x nd x ^ V ‘Zdx^dg^ 

is automatically of first order in d\- Furthermore, we Taylor expand the effective masses 
on the left hand side as well as the propagators and self energies on the right hand side to 
first order in dx around X. After that, we Fourier transform the difference equation with 
respect to s. The Wigner transformed statistical propagator and spectral function are given 
by 

G f ( X , k) — J d 4 s exp (— iks ) G F ( X , s') 

and 

G e (A", k) = —i f d 4 s exp (—iks) G Q (A, s) . 
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Kadanoff-Baym 


Kadanoff-Baym 



FIG. 4: The modulus of the unequal-time FIG. 5: The equal-time propagator as a func- 

propagator as function of time for fixed mo- tion of time for three different initial condi- 

menturn mode p = 0. Correlations be- tions (cf. Fig. 6). The system shows rapid os- 

tween earlier and later times are exponen- dilations which die out after moderate times 

tially damped. and are followed by a smooth drifting regime. 

The factor of — i in the Wigner transform of the spectral function makes G Q (.A, A) again a 
real-valued function. However, in order to be able to really perform the Fourier transfor¬ 
mation, we have to send the initial time to — oo. At least for large x° and y° this can be 
justified by taking into account that correlations between earlier and later times are sup¬ 
pressed exponentially, as one can see in Fig. 4. The result of all these transformations is a 
quantum kinetic equation for the statistical propagator 4 [14-16, 39-42]: 

(2 k»d x , - {d x ,M 2 (A')) d kti ) G f (A, A) 

= fl e (A, k) G f (A', k) - n F (A, k) G e (A, k) (11) 

+ {fl F (A, k ); Re (G R (A, A)) + {Re (fl R (A, A)) ; G F (A, A)}^ , 

where the Poisson brackets are defined by 

{f{X,k)-,g{X,k)} pB = [d x J(X,k)] [d k J (A, A) ] - [d k J(X,k)] [d x „~g(X,k)] . (12) 

Employing the first order Taylor expansion is clearly not justifiable for early times when the 
equal-time propagator is rapidly oscillating, cf. Fig. 5. But this is obvious, since employing 
this gradient expansion is clearly motivated by equilibrium considerations: In equilibrium the 
propagator depends on the relative coordinates only. There is no dependence on the center 
coordinates, and one may hope that there are situations where the propagator depends only 

4 The retarded propagator Gr (x , y) = 9 (x° — y° ) G e (x , y) and self energy, as well as the corresponding 
advanced quantities, have to be introduced in order to remove the upper boundaries of the memory 
integrals. As a result the complete system of quantum kinetic equations includes six equations: one 
equation for G F , G g , Gr , Ilj?, II e and 11#, respectively. 
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moderately on the center coordinates. This is clearly the case for late times when our system 
is sufficiently close to equilibrium. However, as is shown in Fig. 5, already after moderate 
times the rapid oscillations mentioned above, have died out and are followed by a smooth 
drifting regime [17]. In this drifting regime the second derivative with respect to X should 
be negligible as compared to the first order derivative and a consistent Taylor expansion can 
be justified even though the system may still be far from equilibrium. However, it is crucial 
that the Taylor expansion is performed consistently for two reasons: First, this guarantees 
that the quantum kinetic equations satisfy exactly the same conservation laws as the full 
Kadanoff-Baym equations do [15]. Second, it has been shown that neglecting the Poisson 
brackets severely restricts the range of validity of the quantum kinetic equations [42, 43]. For 
the intermediate and late-time regimes these quantum kinetic equations have the advantage 
that they do not include any memory integrals. Being local in time, their numerical solution 
requires much less computer memory as compared to the Kadanoff-Baym equations and 
algorithms using an adaptively controlled time-step size become available. Furthermore, the 
energy convolutions replacing the memory integrals can be done quite efficiently using a Fast 
Fourier Transform algorithm. In order to derive a Boltzmann equation from the quantum 
kinetic equation (11), first we have to discard the Poisson brackets, thereby sacrificing the 
consistency of the gradient expansion. After that, we employ the Kadanoff-Baym ansatz 

G f (X, k ) = G e (X, k ) (n (X, k) + ^j , (13) 

which also can be motivated by equilibrium considerations. In fact, this is a generalization 
of the fluctuation-dissipation theorem, which states that, for a system in thermal equilib¬ 
rium, the statistical propagator is proportional to the spectral function. The fluctuation 
dissipation theorem can be recovered from Eq. (13) by discarding the dependence on the 
center coordinate X and fixing n to be the Bose-Einstein distribution function. The last 
approximation, which is necessary to arrive at the Boltzmann equation, is the so-called 
quasi-particle (or on-shell) approximation: 

G, (X, k) = —£ (s (k° - E (X, k)) - 6 (k° + E (X, fc)) ) , (14) 

where the quasi-particle energy is given by 

E(X,k)=sJm? h + kK 

Once more, we would like to stress that the exact time evolution of the spectral function is 
determined by the Kadanoff-Baym equations. It has been shown that the spectral function 
can be parameterized by a Breit-Wigner function with a non-vanishing width [18, 23]. To 
reduce the width of this Breit-Wigner curve to zero is certainly not a controllable approx¬ 
imation and leads to significant qualitative discrepancies between the results produced by 
Kadanoff-Baym and Boltzmann equations. In fact this approximation can only be justified 
if our system consists of stable, or at least very long-lived, quasi-particles, whose mass is 
much larger than their decay width. We also would like to note that a completely self- 
consistent determination of the thermal mass in the framework of the Boltzmann equation 
requires the solution of an integral equation for E(X,k), which would drastically increase 
the complexity of our numerics. As none of our physical results depend on the exact value of 
the thermal mass, for convenience, we set m t h to the equilibrium value of the thermal mass 


9 



as determined by the Kadanoff-Baym equations. Eventually, we define the quasi-particle 
number density by 

n{X,k) = n{X,k,E{X,k)) . 

After equating the positive energy components in Eq. (11) we arrive at the Boltzmann 
equation. For a spatially homogeneous system there is no dependence on the spatial center 
coordinates and the Boltzmann equation reads 5 : 


d t n (t, k ) = 


A 2 7r 
"48" 


f d 3 p 

/ - q 

f di \ [ <Pr \ 

/ (27r) 3 J 

(27T) 3 J 


Ep ~ Eg 


E k Ep E q E r 
E r ) 


x 6 (k + p — q — r)5 (E k 
x ((1 + n k ) (1 + n p ) n q n r - n k n p (1 + n q ) (1 + n r ) j 


(15) 


For a spatially homogeneous and isotropic system we can dramatically simplify the collision 
integral [44], which allows us to reduce the complexity of our numerics significantly. The 
details of this calculation are shown in the appendix. The result is the following equation 6 : 


d t n ( t , k ) = 


A 2 

967T 4 


OO OO 

I dp ! dq 

0 0 


0 ( r o) 


pqD(k,p,q,r 0 ) 

E k E p Eg 


x ((1 + n k ) (1 + n p ) n q n ro - n k n p (1 + n q ) (1 + n ro ) ) 


(16) 


The auxiliary functions ro and D are obtained from very simple expressions which are given 
in the appendix. In this section we have shown that, using a gradient expansion and a 
quasi-particle (or on-shell) approximation, one can derive the Boltzmann equation from the 
Kadanoff-Baym equations. In this sense one can consider the Kadanoff-Baym equations as 
quantum Boltzmann equations re-summing the gradient expansion up to infinite order and 
including off-shell and memory effects. In the next section we are going to explain how we 
solved the Boltzmann and the Kadanoff-Baym equations numerically. 


V. NUMERICAL IMPLEMENTATION 
A. Kadanoff-Baym Equations 

For the numerical solution of the Kadanoff-Baym equations we follow exactly the lines of 
Refs. [17, 45, 46], i.e. for the spatial coordinates we employ a standard discretization on a 
three-dimensional lattice with lattice spacing a s and N s lattice sites in each direction. Thus, 
the lattice momenta are given by 



5 Here, we use the abbreviations E k = sjir\% h + k 2 and n k = n (t , k). 

6 Now, k = |fc| and n k = n ( t , |fe|) 
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where rij , j G {1,2,3}, enumerates the momentum modes in the j -th dimension. As we 
consider a spatially homogeneous and isotropic system, for given times (x°, y°) we only need 
to store the propagator for momentum modes with > n\ > n 2 > n$ > 0. This saves 
us a factor of 48 in memory usage. The discretization in time leads to a history matrix 
H = {0, a t , 2 a t ,..., (N t — 1) a t } 2 . Here a t is the step size and N t is the number of times in 
each time dimension for which we keep the propagator in memory in order to compute the 
memory integrals. This history cut off can be justified by the exponential damping of the 
unequal-time propagator, cf. Fig. 4. Exploiting the symmetry of the statistical propagator 
with respect to the interchange of its time arguments, we only need to store the values 
of the statistical propagator for x° > y°. In very much the same way we can use the 
respective antisymmetry of the spectral function. This saves us another factor of 2 in 
memory usage. The convolutions arising in the computation of the setting-sun self-energies 
are most efficiently computed using a Fast Fourier Transform algorithm for real-valued even 
functions [47]. 

In order to set the scale for the simulations, we use the renormalized vacuum mass ijir. 
The corresponding bare mass mg is obtained through a perturbative renormalization at 
one-loop order of the self energy (tadpole) [48]. We solved the Kadanoff-Baym equations 
numerically on a lattice with N t = 500, N s = 32, a t mR = 0.06, a s rriR = 0.5 and A = 18. So 
far we performed our simulations on a simple desktop PC with a Pentiumd processor and 2 
GByte RAM. However, we would like to note that the numerics can easily be parallelized. 


B. Boltzmann Equation 

As we saw in the previous paragraph, in order to discretize the Kadanoff-Baym equations 
we can rely on the well-defined scheme offered by standard lattice held theory. Unfortu¬ 
nately, the energy conserving 5 function in Eq. (15) prevents us from using these standard 
lattice techniques for the Boltzmann equation. The reason is the following: When inte¬ 
grating over an arbitrary momentum mode in Eq. (15) one has to look for zeros of the 
argument of the energy conserving 5 function with respect to this particular momentum 
mode. These zeros might well fall between two lattice sites. Hence, computing the colli¬ 
sion integral requires the use of interpolation techniques in order to determine the particle 
number distribution for these in-between lattice sites. These interpolation techniques imply 
a continuity assumption for the particle number distribution which contradicts the strict 
lattice discretization as offered by lattice held theory. Apart from this principal obstacle, 
there is also a practical reason which encourages us to use different discretization schemes for 
both types of equations: The collision integral in Eq. (15) is no convolution. Consequently, 
Fast Fourier Transformation algorithms are not applicable, and its numerical computation 
becomes rather expensive. In order to reduce the complexity of our Boltzmann numerics we 
exploited isotropy, which allowed us to simplify the Boltzmann equation analytically and 
lead us to Eq. (16). In the discretized version of the Boltzmann equation (16) the momenta 
are of the form 

a/T2 

Pn = —Ty re • 

We use the same value for a s as for the Kadanoff-Baym equations. This ensures that the 
largest available momentum is the same as for the Kadanoff-Baym equations. Of course, 
N s need not be the same as for the Kadanoff-Baym equations, which just means that we 
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FIG. 6: Initial particle number densities 

against absolute momenta. Shown are the three 
different initial conditions (IC) discussed in the 
text, for which we numerically solved the Boltz¬ 
mann and the Kadanoff-Baym equations, re¬ 
spectively. All initial conditions correspond to 
the same (conserved) average energy density. 
Above that, the initial conditions IC1 and IC2 
also correspond to the same initial average par¬ 
ticle number density. 


approach the physically relevant infinite volume limit independently for both types of equa¬ 
tions. 

In order to compute the collision integral we proceed as follows: For fixed (k,p,q) we 
determine r 0 (the exact definition of r 0 is given in the appendix), which of course need not 
be one of the discretized momenta given above. The function D ( kpp , q, ?’o) can be evaluated 
for any value of (as one also can see in the appendix). To obtain the particle number 
density for an arbitrary 7*0 we use a cubic spline interpolation [49]. Thus, for given ( k,p,q ) 
the integrand is known to any given accuracy and for given k we can simply sum over p 
and q. In order to advance in time we use a Runge-Kutta-Cash-Karp routine with adaptive 
step-size control [49]. 

In order to set the scale for the simulations, again we use the renormalized vacuum mass 
rriR. Our simulations were done with N s = 500, a s m F = 0.5 and A = 18. 


VI. COMPARING BOLTZMANN VS. KADANOFF-BAYM 


We consider three different initial conditions which correspond to the same average energy 
density. Above that, the initial conditions IC1 and IC2 also correspond to the same initial 
average particle number density. The corresponding initial particle number distributions 
are shown in Fig. 6. These particle number distributions can immediately be fed into 
the numerics for the Boltzmann equation. In order to obtain the initial conditions for 
the Kadanoff-Baym equations, we follow Refs. [17, 31]: The initial values for the spectral 
function are determined from the canonical commutation relations. On the other hand, for 
a given initial particle number density, the initial values for the statistical propagator and 
its derivatives are determined according to: 


G F (X 0 ,y 0 ,p) n 0=0 = 


n(t,p) + \ 


L 

[d x oG F (x°,y°,p)] n 0=0 = 0 , 


J i=0 


[d x od y oG F (x°,y°,p)] xO=yO=0 = 


w(t,p) n(t,p ) + 


t =0 


where the initial effective energy density is given by 


(17) 

(18) 
(19) 
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FIG. 7: These plots show the time evolution of the particle number distributions for two 

different momentum modes and all initial conditions (cf. Fig. 6) as determined by the Boltzmann 
and the Kadanoff-Baym equations, respectively. We see that the Kadanoff-Baym equations respect 
full universality, whereas in the case of the Boltzmann equation only a restricted universality is 
maintained, cf. Fig. 8. 


= 0,p) = \Jm 2 R + p 2 . 

Figs. 7 and 8 show the evolution of the particle number distributions for two momentum 
modes and the corresponding equilibrium particle number distributions, respectively, for all 
initial conditions. In the left plots we can see, that the Kadanoff-Baym equations lead to a 
universal equilibrium particle number density. The left plot in Fig. 7 shows that the particle 
number distributions may evolve quite differently for early times 7 . However, respecting 
universality, for any given momentum mode all distributions approach the same late-time 
value. This plot is supplemented by the left plot in Fig. 8. There, one can see that the 
various particle number densities, after equilibrium has effectively been reached, indeed 
completely agree. Hence, this plot proves that we could have shown plots similar to the 
left one in Fig. 7 for all momentum modes. In particular the predicted temperature, given 
by the inverse slope of the line, is the same for all initial conditions. In contrast to this, 
the right plots reveal that the Boltzmann equation respects only a restricted universality. 
In general, e.g. for the initial conditions IC1 and IC3, for any given momentum mode the 
particle number densities will not approach the same late-time value. For both momentum 
modes shown in Fig. 7 a considerable discrepancy is revealed. However, for the special case 
of the initial conditions IC1 and IC2, which, as mentioned above, correspond to the same 
initial average particle number density, the late-time results do agree 8 . 

The reason for the observed restriction of universality can be extracted from Fig. 9. There 


' As we will see, the steep over-shooting of the particle number distribution leads to a quick kinetic equili¬ 
bration, whereas the rather long tail accounts for chemical equilibration. 

8 In Fig. 8 one can see that in the case of the Boltzmann equation there is only one momentum mode for 
which the late-time values of all particle number densities agree, namely the intersection point of the 
lines. However, we could easily have chosen a fourth initial condition for which the late-time result would 
intersect the lines in Fig. 8 in different points. Then there would not be a single momentum mode for 
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FIG. 8: Here, we plotted the equilibrium particle number distributions obtained for times 
when thermal equilibrium has effectively been reached, against the corresponding thermal energy 
densities. The thermal mass is given by the zero mode of the effective kinetic equilibrium energy 
density as determined by the Kadanoff-Baym equations: rrith = u eq (p = 0). For a given initial 
condition, the temperature is given by the inverse slope of the line and the chemical potential is 
obtained from the y-axis intercept divided by — (3. Supplementing Fig. 7 we observe full (restricted) 
universality in the case of the Kadanoff-Baym (Boltzmann) equations. In particlular, the Kadanoff- 
Baym equations lead to a universal temperature T = 1.68 my, and a universally vanishing chemical 
potential. In contrast to this, the Boltzmann equation gives T = 1.52 rrith and p = 0.18 mth for 
the initial conditions IC1 and IC2, and T = 1.32m^ and p = 0.68 mth for IC3. 


we show the time evolution of the total particle number per volume 

Ntot ( t ) = [ n (t, p ) . 

J (27r) 

In general the Kadanoff-Baym equations conserve the average energy density 9 and global 
charges [24-26]. However, as there is no conserved charge in our theory, the total particle 
number need not be conserved. Indeed, the Kadanoff-Baym equations include off-shell par¬ 
ticle creation and annihilation [18]. Consequently, the total particle number may change, 
and in fact approaches a universal equilibrium value. In contrast to this, due to the quasi- 
particle (or on-shell) approximation particle number changing processes are kinematically 
forbidden in the Boltzmann equation. The Boltzmann equation only includes two-particle 
scattering, which leaves the total particle number constant. Of course, this additional con¬ 
stant of motion severely restricts the evolution of the particle number density. Therefore the 
Boltzmann equation cannot lead to a universal quantum thermal equilibrium. Only initial 
conditions for which the average energy density and the total particle number agree from 
the very beginning, lead to the same equilibrium results. 


which the late-time values of all particle number densities agreed. 

9 Concerning our simulations, of course, this only holds up to numerical errors. We have checked that our 
simulations conserve the average energy density up to a numerical uncertainty of 0.2% for the Kadanoff- 
Baym equations and the Boltzmann equation, respectively. 


14 




Boltzmann 


Kadanoff-Baym 



0.33 


0.32 



O 

Z 0.3 


- IC 1 

- IC 2 
IC 3 


0.29 


n oei _i_i_i_i_i_ i _*_i_i_i_ 

u - 8 0 500 1000 1500 2000 2500 

tm R 


FIG. 9: Time evolution of the total particle number. As expected from Ref. [18], the 
Kadanoff-Baym equations include off-shell particle creation and annihilation. As a result the total 
particle number may change with time. In contrast to this the total particle number is strictly 
conserved in the case of the Boltzmann equation. Concerning our simulations, of course, this 
only holds up to numerical errors (< 0.8%). The quantitative disagreement of the total particle 
numbers in both plots can be attributed to the substantial discrepancies in the discretization 
schemes underlying our Boltzmann and Kadanoff-Baym numerics and are of no relevance for the 
purposes of the present work. 


In a system allowing for creation and annihilation of particles, the chemical potential of 
particles, whose total number is not restricted by any conserved quantity, must vanish in 
thermodynamical equilibrium. The chemical potential predicted by the Kadanoff-Baym and 
Boltzmann equations, respectively, is given by the y-axis intercept, extracted from Fig. 8, 
divided by —j3. Using a ruler the reader might convince himself that the Kadanoff-Baym 
equations indeed lead to a universally vanishing chemical potential. In contrast to this, 
even without a ruler one can see that the Boltzmann equation, in general, will lead to a 
non-vanishing chemical potential 10 . 

In this context, Fig. 10 exhibits further interesting results. In the upper left plot one can 
see that the Kadanoff-Baym equations rapidly wash out our tsunami-type initial condition 
IC3. In both plots on the left hand side the double-dashed-dotted lines correspond to the 
particle number distribution at the same time tuiR = 42.4. Thus, in the lower left plot 
one obtains an approximate straight line already after a relatively short period of time, 
indicating a swift approach to kinetic equilibrium. Subsequently, this straight line is tilted 
until it intersects the origin of our coordinate system (full line), corresponding to a vanishing 
chemical potential . However, this approach to full thermodynamical (including chemical) 
equilibrium takes a considerably longer time [9]. In this way, the left plots reveal two distinct 
time scales: a rather fast kinetic equilibration, and a very slow thermodynamical (including 
chemical) equilibration. These two time scales can also be identified in the left plot of Fig. 7 


10 For the initial conditions considered in this work, the Boltzmann equation predicted even a positive 
chemical potential. However, already on very general grounds, one can deduce that the chemical potential 
of bosons has to be negative! 
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FIG. 10: (Missing) separation of time scales. The particle number distribution is shown 
against the equilibrium energy density at various times for initial condition IC3. 


and in Fig. 5. The over-shooting of the particle number density for early times leads to 
the kinetic equilibration. In fact, the double-dashed-dotted lines correspond to the time, 
when the particle number distribution (equal-time propagator) reaches its maximum value 
in Fig. 7 (Fig. 5). Interestingly, although the initial conditions IC1 and IC2 do not show 
this excessive over-shooting, the corresponding particle number distributions (equal-time 
propagators) approach each other on the same time scale, from which on they show an almost 
identical evolution. The following rather long tail, again indicates that full thermalization 
takes place on much larger time scales. The right plot in Fig. 7 shows that the steep initial 
evolution, which is characteristic for the Kadanoff-Baym equations, is absent in the case of 
the Boltzmann equation 11 and that the Boltzmann equation leads only to a gently inclined 
evolution of the particle number distribution. Accordingly, the plots on the right hand side 
of Fig. 10 show that it takes a considerably longer time for the Boltzmann equation to 
reach kinetic equilibrium. As already mentioned above, in contrast to the Kadanoff-Baym 


11 One might be tempted to conclude that the evolution of the particle number distribution is strictly 
monotonous in the Boltzmann case. However, the small dip for the particle number distribution IC2 in 
Fig. 7 shows that this is not necessarily the case. 
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equations, the Boltzmann equation cannot describe the process of chemical equilibration. 
Consequently, the separation of time scales furnished by the Kadanoff-Baym equations is 
absent in the Boltzmann case. 


VII. CONCLUSIONS 

Starting from the 2PI effective action for a scalar <f> 4 quantum held theory we briefly 
reviewed the derivation of the Kadanoff-Baym equations and the approximations which are 
necessary to eventually arrive at a Boltzmann equation. We solved both, the Kadanoff-Baym 
equations and the Boltzmann equation, numerically for spatially homogeneous and isotropic 
systems in 3+1 dimensions without any further approximations. 

We have shown that the Kadanoff-Baym equations respect universality: For systems with 
equal average energy density the late time behavior coincides independent of the details 
of the initial conditions. In particular, independent of the initial conditions the particle 
number densities, temperatures and thermal masses predicted for times when equilibrium 
has effectively been reached coincide. The chemical potentials also coincide and vanish. 
Furthermore, we observed that thermalization takes place on two different time-scales: a 
rather fast kinetic equilibration, and a very slow thermodynamical (including chemical) 
equilibration. 

In general Kadanoff-Baym and Boltzmann equations conserve total energy as well as 
global charges. In the special case of a real scalar <f> 4 quantum held theory the quasi-particle 
approximation implies that the Boltzmann equation additionally conserves the total particle 
number. This additional constant of motion severely restricts the evolution of the system. As 
a result the Boltzmann equation cannot lead to a universal quantum thermal equilibrium. 
The Boltzmann equation respects only a restricted universality: Only initial conditions 
for which the average energy density and the total particle number agree from the very 
beginning, lead to the same equilibrium results. In particular, the Boltzmann equation 
cannot describe the phenomenon of chemical equilibration and, in general, will lead to a 
non-vanishing chemical potential. Due to the lack of chemical equilibration, the separation 
of time scales, which we observed for the Kadanoff-Baym equations, is absent in the case of 
the Boltzmann equation. 

Some of the approximations that lead from the Kadanoff-Baym equations to the Boltz¬ 
mann equation (namely, the gradient expansion, neglecting the Poisson brackets and the 
Kadanoff-Baym ansatz) are clearly motivated by equilibrium considerations. Taking the 
observed restriction of universality into account, we conclude that one can safely apply the 
Boltzmann equation only to systems which are sufficiently close to equilibrium. Accordingly, 
for a system far from equilibrium the results given by the Boltzmann equation should be 
treated with care. 

In the future it will be important to perform a similar comparison between Boltzmann 
and Kadanoff-Baym equations for Yukawa-type quantum held theories including fermions 
and gauge fields. Also a treatment of Kadanoff-Baym equations on an expanding space-time 
should reveal interesting results. This would finally enable one to attack the problem of 
leptogenesis. Independent of the comparison between Boltzmann and Kadanoff-Baym equa¬ 
tions we are looking forward to learn to which extend a full non-perturbative renormalization 
procedure [34, 50-54] affects the results quantitatively. 
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SIMPLIFYING THE BOLTZMANN EQUATION 

The simplification of the Boltzmann equation [44] relies on the Fourier representation of 
the momentum conservation delta function: 


<r (m) = 


d 3 £ 

( 2^7 


exp . 


Using spherical coordinates, we find 

m £ = m £ ^sinr? m sin cos (y? m — (p^) + cos$ m cost^j . 

Now, we consider just the integration over the solid angle. As we integrate over the complete 
solid angle f^, it does not matter in which direction m is pointing. The result will always 
be the same: 

J dffg exp (— im£) = j exp (— im 0 £) , 

where we can choose m 0 = (0,0, m), such that = 0. Now, we can evaluate the 

integral quite easily: 

f f 4tt 

J dfl^ exp (— im^) = J dfl^ exp (—zm^cost^) = —- sin (m^) . (A.20) 

After we have rewritten Eq. (15) using spherical coordinates and inserted the Fourier repre¬ 
sentation for the momentum conservation delta function, we can use Eq. (A.20) to perform 
the integrations over the solid angles. Here it is crucial first to evaluate the integrals over 
kl q and Q r , and to do the integral over fig at last. We find: 


d t n ( t , k ) = 


A 2 


967T 4 


oo oo oo oo 


dp dq dr d£ 


pqr 


d (Ek + Ep — E q — E r ) 

EkE p Eg 


x ^ 2 sin (^) sin (p£) sin (<?0 sin (^) 


x ( (1 + n k ) (1 + n p ) n q n r - n k n p (1 + n q ) (1 + n r ) j 

There are only two more steps to make in order to arrive at Eq. (16). First, we define the 
auxiliary function D ( k,p , q. r ): 


D(k,p,q,r)= J d£-j ^2 8in(k£) sin (p£) sin (q£) sin (rf) . 
o 
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This can easily be evaluated using a computer algebra program. For k > 0 this is 


D{k,p,q,r ) 


- (\k — p — q — r\ — \k + p — q — r 

16 k V 

— \k — p+q — r\ + \k + p+ q — r\ 

— \k — p — q + r\ + \k + p — q + r\ 

+ \k — p + q + r\ — \k + p + q + r\} 


and for k = 0 we obtain 

D(0,p,q,r) 


- ^sign (p + q — r) — sign {p — q — r) 

+ sign (p — q + r) — sign (p + q + r) j. 


Second, we use the energy conservation <5 function to evaluate the integral over r, using the 
well-known formula 




E 

Fol/ho)=o} 




ro) 


(4L) 

\dr ) r=ro 


r 0 is determined by the condition that the argument of the energy conservation 5 function 
is zero: 


Ek + E p — E q — E ro — 0 . 


(A.21) 


If this condition can be satisfied, r' 0 is given by 


r 0 = r 0 (t, k, p, q) = yj(E k + E p - E q ) 1 2 3 4 5 6 7 - M 2 (t) . 

If k, p and q are such that condition (A.21) cannot be satisfied, the above square root yields a 
purely imaginary result and Tq < 0. Due to the 0 function the corresponding term does not 
contribute to the collision integral. After these final steps we end up exactly with Eq. (16). 
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